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Abstract: We revisit the Higgs portal vector dark matter model including a hidden 
sector Higgs field that generates the mass of the vector dark matter. The model becomes 
renormalizable and has two scalar bosons, the mixtures of the SM Higgs and the hidden 
sector Higgs bosons. The strong bound from direct detection such as XENONIOO is evaded 
due to the cancellation mechanism between the contributions from two scalar bosons. As 
a result, the model becomes still viable in large range of dark matter mass, contrary to 
some claims in the literature. The Higgs property is also affected, the signal strengths for 
the Higgs boson search being universally suppressed relative to the SM value. 
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1 Introduction 

The so-called Higgs portal cold dark matter (CDM) model is an interesting possibility for 
the nonbaryonic dark matter of the universe. The dark matter fields are assumed to be 
the SM gauge singlets, and could be a scalar (S), a singlet fermion (ip) or a vector boson 
(y) depending on their spin. The Lagrangian of these CDM's are usually taken as [1] 

/:.caiar = Id.Sd'^S - ImlS' - ^H^HS' - ^S' (1.1) 

Aermion = 1p [i^f- 8 - 171^] l/j - -j^H^H iplp (1.2) 

/^vector = V^''' + l^vy^^y^ + |Ay(F^y^)2 + ^Xhv HV^V^' . (1.3) 

Dark matter fields {S, ip, V) are assumed to be odd under new discrete Z2 symmetry: 
(5, ■0, y) —{S, ip, V) in order to guarantee the stability of CDM. This symmetry removes 
the kinetic mixing between the V^u and the U{1)y gauge field B'^'^ , making V stable. 

The scalar CDM model (1.1) is fine, as long as Z2 symmetry is unbroken. The model 
is renormalizable and can be considered to high energy scale as long as the Landau pole is 
not hit. On the other hand, the other two cases have problems. 

Let us first consider the fermionic CDM model (1.2). This model is nonrenormalizable, 
and has to be UV completed. The simplest way to achieve the UV completion of (1.2) is to 
introduce a real singlet scalar field as proposed in Ref. [2, 3] by some of us. We observed that 
there are two Higgs- like scalar bosons which interfere destructively in the spin-independent 
cross section of the singlet fermion CDM on nucleon. The strong constraint from direct 
detection experiments such as XENONIOO [4] or CDMS [5] can be relaxed by a significant 
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amount, unlike the claim made in literatures [1] based on the effective Lagrangian (1.2). 
The decoupling of the 2nd scalar boson occurs rather slowly, since the mass mixing between 
the SM Higgs boson and the new singlet scalar is due to the dim-2 operator. Also the mixing 
between two scalar bosons makes the signal strength of two physical Higgs-like bosons less 
than one, and make it difficult to detect both of them at the LHC. Since there is now an 
evidence for a new boson at 125 GeV at the LHC [6, 7], the 2nd scalar boson in the singlet 
fermion DM model is very difficult to observe at the LHC because its signal strength is 
less than 0.3 [3, 8]. Also an extra singlet scalar saves the vacuum instability for rrin = 125 
GeV [8-10]. The electroweak (EW) vacuum can be still stable upto Planck scale even for 
niH = 125 GeV [8]. These phenomena would be very generic in general hidden sector DM 
models [11]. Li short, it is very important to consider a renormalizable model when one 
considers the phenomenology of a singlet fermion CDM. 

Now let us turn to the Higgs portal vector dark matter described by (1.3) [1]. This 
model is very simple, compact and seemingly renormalizable since it has only dim-2 and 
dim-4 operators. However, it is not really renormalizable and violates unitarity, just like the 
intermediate vector boson model for massive weak gauge bosons before Higgs mechanism 
was developed. The Higgs portal VDM model based on (1.3) is a sort of an effective 
lagrangian which has to be UV completed. It lacks including the dark Higgs field, (^(x), 
that would mix with the SM Higgs field, h{x). Therefore the model (1.3) does not capture 
dark matter or Higgs boson phenomenology correctly. It is the purpose of this work to 
propose a simple UV completion of the model (1.3), and deduce the correct phenomenology 
of vector CDM and two Higgs-like scalar bosons. Qualitative aspects of our model are 
similar to those presented in Ref.s [3, 8], although there are some quantitative differences 
due to the vector nature of the CDM. 

This work is organized as follows. In Sec. 2, we define the model by including the 
hidden sector Higgs field that generates the vector dark matter mass by the usual Higgs 
mechanism. Then we present dark matter and collider phenomenology in the following 
section. The vacuum structure and the vacuum stability issues are discussed in Sec. 4, and 
the results are summarized in Sec. 5. 

2 The model 

Let us consider a vector boson dark matter, X^, which is assumed to be a gauge boson 
associated with Abelian dark gauge symmetry U{l)x- The simplest model will be without 
any matter fields charged under U{l)x except for a complex scalar, $, whose VEV will 
generate the mass for [12]: 



^VDAI 




(2.1) 



in addition to the SM lagrangian. The covariant derivative is defined as 



D^<^ = {d^ + igxQ^X^)<^, 
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where Q$ = Qx{^) is the U{l)x charge of ^ and we wih take Q$ = 1 throughout the 
paper. 

Assuming that the C/(l)x-charged complex scalar $ develops a nonzero VEV, and 
thus breaks U{l)x spontaneously, 

^ = ^ (^^K + ^i^)) ■ 

Therefore the Abelian vector boson get mass = gxQ%^^^ and the hidden sector 
Higgs field (or dark Higgs field) ^p{x) will mix with the SM Higgs field h{x) through 
Higgs portal of the Xh^ term. It is straightforward to derive the Feynman rules from 
this Lagrangian after the U{l)x and EW symmetries are spontaneously broken. 



3 Phenomenology 

3.1 Dark Matter Phenomenology 

The observed present cold dark matter density, UcDuh'^ — 0.1123 it 0.0035 [13], is approx- 
imately related to the thermally averaged annihilation cross section at freeze-out temper- 
ature, {crv){z, as 

2 _ 3 X 10-^W/s 

So we require (<Tf)fz 3 x 10~^^cm'^/s to obtain the correct relic density. We have used 
the micrOmegas v. 2. 4. 5 [14] to calculate thermal relic density and direct detection cross 
section of the VDM in our model. 

In Fig. 1 we show the thermal relic density as a function of the dark matter mass, Mx- 
For this plot we fixed mi = 125 GeV , m2 = 150 GeV , a = 7r/4 and the purple (blue) line 
corresponds to gx = 0.05 (0.5). We can see two resonance dips at Mx = rrii/l {i = 1,2). 
The VDMs can annihilate into the SM particles in the S-wave state, which is different 
from the singlet fermionic dark matter case studied in [3] where the annihilation occurs in 
the P-wave state. As a result the annihilation cross section for the vector dark matter is 
generally 0(10 — 100) larger than that of the SFDM. And the current relic density can be 
explained more easily even at non-resonance region. (See the blue line in Fig. 1.) 

One important effect when considering the full theory, which we found in Ref. [3], is 
that a generic cancellation occurs in the dark matter and nucleon scattering amplitude, 
which can not be observed in the effective lagrangian approach. ^ This is because the 
transformation matrix between the interaction eigenstates and the mass eigenstates in the 
scalar sector is an orthogonal matrix. The dark matter and nucleon elastic scattering cross 

^In general the cancellation mechanism can also work in the annihilation process for the relic density. 
However, the difTerent decay widths for the Hi and H2 and/or other processes such as annihilations into 
scalar particle pairs makes it less effective than in the direct detection process. As a result, the annihilation 
process and the direct detection process are not strictly proportional to each other in our scenario. 
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Figure 1. The thermal rehc density fix^i^ of the vector dark matter as of function of the dark 
matter mass, Mx- For this plot we fixed mi = 125 GeV , m2 = 150 GeV ,a = 7r/4 and the purple 
(blue) line corresponds to gx — 0.05 (0.5). The horizontal line is the central value of the current 
relic density ^xh'^ = 0.1123 [13]. 



section is proportional to the following factor^: 



i=l,2 ^ 



OhiO. 



(3.3) 



where q is the momentum transfer of the dark matter. When nii ~ m2 or \q^\ ^ mf, 
we have dp ~ due to the orthogonality of the mixing matrix O. This cancellation phe- 
nomenon is quite similar to the GIM-mechanism [15] in the quark (or lepton) flavor violating 
neutral current processes. In Fig. 2, we show the excluded region in the (g'x, a)-plane by 
the non-observation of dark matter by the XENONIOO which currently gives the strongest 
bound on the dark matter direct detection cross section [4] . Each colored region is excluded 
by XENONIOO direct detection experiment for the m2 value given in the plot. We fixed 
Mx = 70 GeV, mi = 125 GeV for the plot. The black line corresponds to 17x/i^ = 0.1123 
with m2 = 135 GeV. So the relic density is consistent with the XENONIOO experiment. 
The entire region is also allowed by the electroweak precision 5, T, [/-parameters at 99% 
confidence level for this choice of parameter space. 

The predictions of our model on the S", T parameters assuming [/ = are shown 
in Fig. 3 for the choices (mi, ma) = (25,125), (50,125), (75,125), (100,125), (125,125), 



^The transformation matrix O is defined as 



= O 



Hi 



Hi 

H2 



(3.2) 



where Sa{ca) = sinQ(cosa?), h,ip are the interaction eigenstates and Hi{i = 1,2) are the mass eigenstates. 
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Figure 2. The excluded region in the {gx, a)-plane. Each colored region is excluded by XENONIOO 
direct detection experiment for the TO2 value given in the plot. We fixed Mx = 70GeV,mi = 
125 GeV. The black line corresponds to ^Ixh^ = 0.1123 with 1112 — 135 GeV. So the relic density 
is consistent with the XENONIOO experiment. The entire region is also allowed by the S,T,U- 
parameters at 99% confidence level. 



(125,250), (125,500), (125,750) (GeV). The green (red) dots are for a = 45°(20°). The 
thick black line is the prediction of the SM with the mn in the range [125, 720] (GeV). The 
ellipses represent 68, 95, 99% CL experimental lines from inside out. 

3.2 Collider Phenomenology 

Since the scalar sector is extended, the Higgs phenomenology is different from that of the 
SM. In this subsection we study the possibility that the second Higgs which our model 
predicts could be discovered at the LHC. We will also see that the combination of the 
collider signatures and the DM direct searches is robust enough to exclude or confirm our 
model in the on-going LHC and the next generation DM direct detection experiments. 
The signal strength of a scalar boson Hi=i^2 defined as 

_ ajpp ^ H,)B{H, ^ X) 

can be measured at the LHC. Here i = 1,2 and X is a specific SM final state which the 
scalar boson Hi can decay into. In our model it can be written in terms of T*°*'^^ (^ = 1^ 2) 
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Figure 3. The predictions of (5, T)-paranieters in our model for {mi, 1712) — (25,125), (50,125), 
(75, 125), (100, 125), (125, 125), (125, 250), (125, 500), (125, 750) (GeV). The green (red) dots are for 
a = 45° (20°). The thick black line is the prediction of the SM with the tuh in the range [125, 720] 
(GeV). The ellipses represent 68, 95, 99% CL experimental lines from inside out. 

which is the total decay width of Hi in the SM assuming Hi is a pure SM Higgs and r*°* 
which is the total decay with of Hi in our model [3, 8]: 

ptot.SM 

r. = Ot^—^, (3.5) 

i 

where Ohi = Ca, Oh2 = Sa- The total decay widths can be decomposed as 

ptot _ 2 -ptotjSM . 2 -ptot.hid 
-■-1 — C„i ]^ -I- S^l I , 

rt°* = r*°*'^^ + clT'°'''^'^ + T{H2^ HiHi), (3.6) 

where p*°*''^^'^ ig the total decay width of Hi into the hidden sector assuming Hi is a pure 
SM-singlet scalar. The channel H2 — s- HiHi opens when m2 > 2mi. From the eqs. (3.5) 
and (3.6) it is obvious that rj < 1 in our model. Therefore if the excess of the signal 
strength in some channels like — )• 77 above the SM prediction at the LHC remains 
in the future data, our model will either be excluded or need to be extended (two Higgs 
doublet portal to a hidden sector dark matter, for example). In fact we find r < 0.3 for 
the second Higgs boson, when we impose the observed new boson at 125 GeV as one of 
the Higgs bosons in our model. 
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Figure 4. The scattered plot in (a) (ri,r2) for mi(= 125 GeV) < m2 and (b) (r2,ri) for mi < 
m2(= 125 GeV). The big (small) points (do not) satisfy the WMAP relic density constraint within 
3 (7, while the red- (blue-) colored points can (cannot) be probed at the planned XENONIT direct 
detection experiment. 



The correlation between ri and r2 can be seen in Fig. 4 where we show only the region 
ri > 0.7. For this plot we scanned the parameters gx, Mx, a, fn2 in the range, Q < gx < 
10 GeV < Mx < lOOOGeV, -7r/2 < a < ir/2, mi{= 125 GeV) < m2 < 2000 GeV for the 
panel (a), and 10 GeV < mi < m2{= 125 GeV) for the panel (b). All the points pass the 
constraints: Qxh"^ < 0.1228 (the 3a upper bound of the relic density), the upper bound on 
the XENONIOO direct detection cross section, and the bound on the S, T-parameters at 
99% CL. The big (small) points (do not) satisfy the WMAP relic density constraint within 
3 a, while the red-(blue-)colored points can (cannot) be probed at the planned XENONIT 
direct detection experiment [16]. In both plots, the big red points on the straight line. 
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The 2nd Higgs mass(GeV) 

Figure 5. The allowed mixing angle a as a function of the second Higgs mass. We fixed the 
SM-like Higgs mass to be 125 GeV. The big (small) points (do not) satisfy the WMAP relic density 
constraint within 3cr, while the red-(blue-)colored points can (cannot) be probed at the planned 
XENONIT experiment. The black points are excluded by the LHC Higgs search, i.e. r < 0.7. 

ri + r2 = 1, are those with Hi — XX and H2 — ?■ HiHi suppressed. In the panel (a), 
the sizable contribution from the H2 — >• HiHi channel allows the big red points below the 
ri + r2 = 1 line. 

In Fig. 5, we show the allowed mixing angle a as a function of the second Higgs mass. 
We fixed the SM-like Higgs mass to be 125 GeV. The big (small) points (do not) satisfy the 
WMAP relic density constraint within 3a, while the red- (blue-) colored points can (cannot) 
be probed at the planned XENONIT experiment. The black points are excluded by the 
LHC Higgs search, i.e. r < 0.7. We can see the maximal mixing angle a = 7r/4 (black 
points near 7712 ~ 125 GeV) is excluded by the LHC Higgs search. Also the light scalar 
with mass less than 125 GeV, if exists, should be singlet-like. 

In Fig. 6, we show a scattered plot of (Jp 3(S Si function of Mx- The big (small) 
points (do not) satisfy the WMAP relic density constraint within 3 a, while the red- 
(black-)colored points gives ri > 0.7(ri < 0.7). The Gray region is excluded by the 
XENONIOO experiment. The dashed line denotes the sensitivity of the next XENON 
experiment, XENONIT. We note that many points are still allowed by the WMAP relic 
density constraint, the XENONIOO direct detection experiment, and also by the constraint 
ri > 0.7 which is in the ball park of the LHC Higgs search bound. On the other hand, the 
effective field theory approach considered in Ref.s [1] strongly constrains the vector dark 
matter scenario. We can also see that there is no point below about Mx ~ 50 GeV in 
Fig. 6 (a). It is because the Higgs exchanged dark matter annihilation channel does not 
allow the resonance and the relic density is larger than the WMAP measurement. Most of 
the big red points are within the reach of the XENONIT sensitivity, and our model can 
be tested in the next generation dark matter detection experiment. 
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Figure 6. The scattered plot of dp as a function of Mx- The big (smaU) points (do not) satisfy the 
WMAP rehc density constraint within 3 tr, while the red-(black-)colored points gives ri > 0.7(ri < 
0.7). The gray region is excluded by the XENONIOO experiment. The dashed line denotes the 
sensitivity of the next XENON experiment, XENONIT. 

4 Vacuum stability and perturbativity of Higgs quartic couplings 

In this section, we analyze vacuum stability and perturbativity of Higgs quartic couplings. 
To make the Higgs potential be bounded-from-below, we require 

Xh>0, A$>0, -2^/\hX<s, < \h<s>, (4.1) 
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where the last condition apphes for Xh<^ < 0. We also require 



det M^^^g, = det ^^^f^^^ ^^0^* j = (4AhA^ - A?,^)44 > 0. (4.2) 

Since there is additional direction of the Higgs potential can have minima other than 
our EW vacuum. In the following, we investigate whether the EW vacuum is global or not. 
We closely follow the analysis done in Ref. [8]. 

The tree-level effective potential takes the U{\)x symmetric form 

(4.3) 

where ^pH and ip^ are spacetime-independent classical fields. Unlike the general Higgs 
potential, only nontrivial phase may be the I~phase which is defined by a minimum of 
Vb(0, f^). Such a minimum is given by 



= Mhl + ^^1,. (4.4) 



The differences of vacuum energies of the I- and the EW phases is 

T.(I)/n - x T/(EW)/ ^ 4 , ^H'S> 2 2 A$ 4 4 

^ {4.XHX<s>-Xjj^)vjj, (4.5) 



16A$ 

where we have used Eq. (4.4) in the second line. Therefore, as long as Eqs. (4.1) and (4.2) 
are satisfied, the EW vacuum is always the global minimum. Note that this is not the case 
for the generic Higgs potential [8]. 

Although the EW vacuum is stable at the EW scale, its stability up to Planck scale 
(Afpi ~ 1.22 X 10^^ GeV) is nontrivial question since a renormalization group (RG) effect of 
the top quark can drive Xh negative at certain high-energy scale, leading to an unbounded- 
from-below Higgs potential or a minimum that may be deeper than the EW vacuum. We 
will work out this question by solving RG equations with respect to the Higgs quartic 
couplings and the U{l)x gauge coupling. The one- loop f3 functions of those couplings are 
listed in Appendix A. In addition to the vacuum stability, we also take account of the 
perturbativity of the couplings. To be specific, we impose Xi{Q) < An {i = H, H^, <I>) and 
gj^iQ) < 47r up to Q = Mpi. 

Fig. 7 shows the vacuum stability and the perturbativity constraints in the a-m2 plane. 
We take nii = 125 GeV, gx = 0.05, Mx = m2/2 and v<j, = Mx/{gxQ<^)- The vacuum 
stability constraint is denoted by red line; i.e., the region above the red line is allowed 
for a > 0, and it is the other way around for a < 0. The perturbativity requirement is 
represented by blue line; i.e., the region below the blue line is allowed for a > 0, and it is the 
other way around for q < 0. For a < 0, the region above the dotted black line is excluded 
by Eq. (4.1). Putting all together, for a > the region between the red and blue lines 
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Figure 7. The vacuum stability and perturbativity constraints in the a-m2 plane. We take 
TOi = 125 GeV, gx = 0.05, Mx = m2/2 and v$ = Mx/{gxQ<^)- 

is allowed while for a < the region between the dotted black and blue lines is allowed. 
It should be noted that since the coefficient of A//$ in is doubled in comparison with 
the real singlet case, the improvement of the vacuum stability by the increase of Ajy$ or, 
equivalently a, is more effective. However, unlike the general Higgs potential involving 
explicit U{l)x breaking terms, the EW vacuum cannot be stable up to Planck scale if a 
is exactly zero. 

In Fig. 8, we show the vacuum stability and perturbativity constraints in the Mx-'m2 
plane. We fix a = 0.1 varying g-^, i.e., g-^ = 0.1 (Left Panel) and 0.5 (Right Panel). Once 
gj^ is fixed, the small Mx is realized by a small In such a case, the large m2 is possible 
only by a large A$ since m2 — \/2X^v^ for a small a. This explains the regions excluded 
by A$(Q) > Att in both plots. Indeed, the g^ = 0.5 case yields the severer constraints. As 
for the vacuum stability constraint, the change of g^ has little effect on it, which can be 
understood from the expression of /3ah' Eq. (A. 2). 

5 Conclusions 

In this paper, we revisited the Higgs portal vector dark matter including the hidden sector 
Higgs field $ that provides the vector dark matter mass. Including the hidden sector Higgs 
field makes the model renormalizable and unitary. The constraint from direct detection 
cross section (XENONIOO) still allows a large parameter space in this model. On the 
contrary to some claims that the Higgs portal dark matter model is strongly constrained 
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Figure 8. The vacuum stability and perturbativity constraints in the Mx-m2 plane. We set 
— 0.1 (Left Panel) and 0.5 (Right Panel) with being a = 0.1. 



by XENONIOO data, we showed that the model is still viable. It is crucial to work with 
a model that is renormalizable, and not with effective lagrangian, as in the Higgs portal 
fermion DM model in Ref. [3, 8] Including the hidden sector Higgs field also improves the 
vacuum stability of the model for mn = 125 GeV upto the Planck scale as in Ref. [8]. Our 
model can be tested at colliders by searching for the 2nd Higgs boson and/or the signal 
strength of the 125 GeV Higgs boson. It would take long in order to observe the 2nd 
Higgs boson since its signal strength is smaller than 0.3. In our model, is universally 
suppressed relative to the SM case for all channels. This could be a useful criterion when 
the signal strengths of 125 GeV Higgs boson are measured with smaller uncertainties. If 
Tj is not universally suppressed or larger than one, then our model shall be excluded. 
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A One-loop /3 functions of Higgs quartic couplings 

The renormalization group equation and the f3 functions are given by 

d\{t) 



dlog(Q) 



/3a, (A.l) 
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where 



1 



167r2 
1 

167r2 
1 



+ - 6yf + ^{2gl + {gl + 5?)^} - XH{3{3gl + 5?) - 12y?}l , 



(A.2) 
, (A.3) 



167r2 

1 1 3^2 



2Ah$(6Aj^ + 4A$ + 2Xh^) - Xh^ !^^{3gl + gj) - + Ggj^Ql 
2{Xl^ + 10A| + 3gj,Qi) - UX^g^Ql] , (A.4) 



Pgx = i^^aScQi- (A.5) 
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